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Abstract — In this work we consider the problem of 
downlink resource allocation for proportional fairness 
of long term received rates of data users and quality 
of service for real time sessions in an OFDMA-based 
wireless system. The base station allocates available 
power and bandwidth to individual users based on 
long term average received rates, QoS based rate 
constraints and channel conditions. We solve the 
underlying constrained optimization problem and pro- 
pose an algorithm that achieves the optimal allocation. 
Numerical evaluation results show that the proposed 
algorithm provides better QoS to voice and video 
sessions while providing more and fair rates to data 
users in comparison with existing schemes. 

I. Introduction 

Recent IEEE 802.16 standard defines the air in- 
terface and medium access control (MAC) speci- 
fications for wireless metropolitan area networks. 
Such networks intend to provide high speed voice, 
data and on demand video streaming services for 
end users. IEEE 802.16 standard is often referred 
to as WiMax and it provides substantially higher 
rates than cellular networks. Besides it eliminates the 
costly infrastructure to deploy cables, therefore it is 
becoming an alternative to cabled networks, such as 
fiber optic and DSL systems [1],[2]. 

Transmissions in IEEE 802.16 networks is based 
on Orthogonal frequency division multiplexing 
(OFDM), which is a multicarrier transmission tech- 
nique that is proposed for high speed wireless trans- 
mission. It is based on a large number of orthogonal 
subchannels, each working at a different frequency. 
Within OFDMA framework, the resource allocated 
to the users come in three dimensions: time, fre- 
quency and power. This requires the scheduler to 
operate with higher degree of freedom but also 
makes the notion of resource fairness obsolete and 
makes the problem more involved. 



There are three main issues that need to be 
considered in multiple access resource allocation, 
which are spectral efficiency, fairness and Quality of 
Service (QoS). These are often contradicting aims 
and achieving these in an OFDMA based system 
is a challenging task. OFDMA-based scheduling for 
systems with heterogeneous traffic has not been stud- 
ied sufficiently. Recently in [3], [4], [5], proportional 
fair scheduling was studied for multicarrier systems. 
However, the scheduling rules in these works do 
not apply sufficiently to different QoS requirements 
and heterogeneous traffic. Besides they only provide 
fairness for short term received rate, instead of long 
term. In this work our goal is to find algorithms 
that provide proportional fairness to data users and 
QoS for real time sessions such as voice and video 
streaming. 

II. System Model 

We consider a downlink system, where a base 
station transmits to N users. Each user is assumed 
to demand a single type of traffic. Let W and P 
denote the total bandwidth and power, respectively. 
Total bandwidth W is divided into N su b subchannels 
of length W su b, each consisting of a group of carriers. 
We consider PUSC as the mode of subchanneliza- 
tion, which is more suitable for mobile users than 
AMC It provides frequency diversity and inter-cell 
interference averaging. This minimizes the perfor- 
mance degradation due to fast fading characteristics 
of mobile environments. We assume flat fading, 
which is a reasonable assumption with PUSC. 

The noise and interference power density is No, 
and the channel gain averaged over the entire band 
from the BS to user i at time t is hi(t), where /?,(/) 
includes path loss, shadowing (lognormal fading) 
and fast fading. There are three classes of users. 
Users in the classes Uo, Us and Uy demand data, 



voice and video traffic, respectively. The system that 
we consider is time slotted with frame length T s . 
The scheduler makes a resource allocation decision 
at each frame. Active period in a voice conversation, 
streaming duration and file size are both very long 
with respect to the frame size. Therefore it is realistic 
that during the course of optimization the number of 
active sessions are fixed. 

In problem formulation we will adopt the Shannon 
channel capacity for AWGN channel as a function 
for bandwidth and transmission power assigned to 
user i: 

,Pi(t)hi(t) 



r i (wi(t),p i (t))=Wi(t)log[ 1 + p 



(1) 



N Wi(t) 

The reason for using Shannon capacity is its simplic- 
ity, and it also approximates the above set of rate- 
SINR relation with P = 0.25. The parameter < 
P < 1 compensates the rate gap between Shannon 
capacity and rate achieved by practical modulation 
and coding pairs listed in IEEE 802.16a/e standards. 
After finding power and bandwidth pairs by solving 
the convex programming problem, we quantize the 
SINRs to the thresholds determined in the standard. 
In this paper, we propose a Fair and QoS-based 
Power and Subchannel Allocation (FQPSA) Algo- 
rithm. Objective of the algorithm is to maximize 
proportional fairness for data users subject to rate 
constraints of real time users selected each frame 
according to a metric. 

III. Problem Formulation 

A. Data Traffic 

The objective for the data users is to optimize 
log sum of the exponential averaged rates of the 
users Average received user rates are updated at 
each frame according to the exponential averaging 
formula: 

Ri(t) = afii (t - 1 ) + ( 1 - Oj) n ( Wi (t ) , Pi (t) ) , V/, t (2) 

This way we consider both current rate as well as 
rates given to the user in the past. Observed at time t, 
the highest consideration is given to the current rate 
r(t), and the rates received at the past t—\,t — 2,... 
carry diminishing importance. 

C(R(0) = flog^(0 = 

i=i 

C(R(>-l)) + £lc(a + Br '^y» ) (3, 



As a matter of fact, we limit ourself to greedy 
schemes in the sense that at frame t, we try to 
maximize the proportional fair capacity C(R(t)) 
without considering the future frames t + lj + 2, 
etc. Only the second term in (f3]) needs consideration. 
Here the log-sum can be replaced with product. 

B. Real Time Traffic 

The Base Station first chooses a number of voice 
and video users according to a metric that reflects 
current channel condition, head of line packet delay 
and long term received rate. Then for the selected 
users the number of bits to be transmitted in the 
current frame is determined. These procedures are 
explained in [6], in detail. Let U' R be the set of 
real time sessions selected and let r\ be the rate 
requirement of real time session i. 

The resulting optimization problem considering 
all types of sessions is as follows: Find, Q: 



(p*,w*) = argmax Y[ « 
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P *,w* > o,Vieu D uu R (7) 



where n-, = 



IV. Solution of the Optimization Problem 

The objective function (01) is an increasing func- 
tion of (w,p), therefore the maximum is achieved 
only when the constraints (|4l [5J [6]) are all met 
with equality. For this reason we can replace these 
inequalities with equalities in the discussion below. 

Actually there is no guarantee that a solution 
can be found that satisfies both rate and power 
constraints, given the rate requirements and channel 
conditions. The rate requirement can be too high 
that it may be impossible to satisfy with the given 
channel conditions. To start with, we assume that 
the problem is feasible. We will discuss about how 
to detect infeasibility of the problem and what to do 

'Here pi,wi, «; are the values at time t. The time index is not 
shown for convenience. 



in that case in the next section. We can write the 
Lagrangian of the problem as: 



L(w,p,X p ,X w X)=Y[ ( a ' + 

ieu D 



a,w,log(l + ^- 
Ri 



+K\ p - L Pi)+K[w- £ Wi 

\ ieU D uU' R J \ ieU D uU' R , 



+ w/log 1 + 
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From now on we assume that all rates are in nats/sec 
and logarithms are natural. Taking the derivatives of 
L(p,v/,X p ,X w ,X r ), we get the following: 

A. For users i £ t/a 



dL(p,v/,X p ,X w X) 



dpi 
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where a, = -j-^-. Let A p = |- and A vl . = B 
dividing (O with ( fTOl we can write for all i G Ud'- 



^ = A fl = n, ((1 +x*)log (1 +**) -jc?) , (1 1) 

where x* = -A? denotes the optimal effective SINR, 
which is the SINR multiplied by the SINR gap 
parameter p\ 

B. For users i 6 t/^ 
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With ^ = ^, dividing equation (O to (O for 
all i £Ujf again gives (fTTb . 

By dividing (fT3l to (fT2l) we can eliminate ?ij"s 
from the problem. It is worth noting that we get the 
same relation between A a /n,- and x, for both data 
and real time sessions. At this point it is convenient 
to define the function f„(x) as: 

/ B (jc) = (l+Jc)log(l+Jc)-Jc. (14) 
Then we have, 

Xi = f; l (K/ni),Vi eu D uu' R (15) 

Signal to noise ratio (jc,-) is a monotonic increasing 
function of A a for all users i € UnUUg. 
For real time users we also have: 



dL(p,yr,X p ,K,V] 
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C. For all nodes i £ Ud U U' r 



dL(p,-w,X p ,X w ,X r 



dX p 

dL(p,w,X p ,X w ,X r ) 
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From Equation (O for data users we can write: 

[A p -n i (l+x*)R i a i } + 
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The 



log(l+*?)(l+**) 
] + operator in Equations (fT9l),(l20l) guarantees 



that Wj,pi > for all users. Given A p and A a we can 
compute the power and bandwidth for users i £ Lf D 
using ( fT9l ) and (l20l) . Given A fl , we can calculate the 
power and bandwidth for users i EU R using ( fT6l ). 
Please note that just from (O, (O and (|20]>, the 
bandwidth and power constraints d5) (0]) are not 
necessarily satisfied. We need to find the right A a 
and A p so that the power and bandwidth constraints 
are satisfied. Let S p (A a ,A p ) and S w (A a ,A p ) be the 
total bandwidth and total power corresponding to A a 
and A p : 



S w (A a ,A p ) = £ Wj(A a ,A p ) 

ieu D uUg 

y [A p -n i (l+x i )R i a i \ + | y 
kj D log{\ +Xi){\+Xi)ni ^logil+xi 
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(21) 



S p (A a ,A p ) = £ pi(A a ,A p ) = 

ieu D yju' R 

y [Ap -nj(\ + x i )R i a i ] + x i y r'jXjjii 
k D log(l +*,)(! +**) j^log(l+x,) 

(22) 

where x-, = f~ l (A a /«,) is the SNR of user i. As a 
result, the problem is finding A* and A* such that 

s w (a;,a;) = w (23) 
s p (K>K) = p (24) 

using Equations £[5]>,(l23]) and ([24]). Note that al- 
though A a and A p are independent variables that 
determine power and bandwidth for each node, they 
become dependent when the power and bandwidth 
constraints d23l (1241 need to be satisfied. 



D. Feasibility of the Solution 

From (fT9l) and (l20l) . A p = corresponds to the 
case that no bandwidth and power is allocated to 
data sessions. If the problem is feasible, then there 
exists a A a that satisfies the following: 

S w (A a ,0) = W, (25) 
S P (A a ,0) < P. (26) 

Let A° be the smallest A a that satisfies (1251 ) and 
(l26l ). In the report [6] we proved the existence 
and uniqueness of A" for a feasible problem using 
monotonicity properties of the functions ^(A^A^) 
and S p (A a ,A p ). 

V. Proposed Algorithm 

Using (fT5l).(l2~Tl) and (l22l) we can develop an 
algorithm to determine the power and bandwidth 
allocation. The algorithm is also able to detect 
infeasibility if there no solution exists. 

Algorithm: 

1) Find A°. If S p (A° a ,0) < P then the problem is 
feasible. If not, find the real time session that 
consumes the largest power and decrease its 
rate. Repeat Step 1 and 2 until feasibility is 
achieved. 

2) Using two nested binary searches find the 
unique (A*,A£'(A*)) that satisfies power and 
bandwidth constraints. 

3) Compute optimal power and bandwidth using 
the optimal (A:,A;(A*)). 

Let us summarize the nested binary search. We 
start from a A a , then we find the A p (A a ) that satisfies 
the bandwidth constraint. If the power constraint 



is exceeded then we decrease A' a , else increase 
it, and search for A*(A^). We continue until both 
constraints are satisfied. 

Since the problem is convex (proved in [6]), 
the problem has a unique solution. We need the 
monotonicity property of function S p (A a , A*(A a )) 
(sum of powers) as a function of A a in order to be 
able to find it using binary search. Those properties 
are also verified in [6]. 

VI. Performance Evaluation 

For the numerical evaluations we equally divide 
the users to 5 classes according to the distances, 
0.3,0.6,0.9,1.2,1.5 km. We use the parameters in 
Table H 



Parameter 


Value 


Cell radius 


1.5km 


User Distances 


0.3,0.6,0.9,1.2,1.5 km 


Total power (P) 


20 W 


Total bandwidth (W) 


8.3 MHz 


Frame Length 


1 msec 


Voice Traffic 


CBR 32kbps 


Video Traffic 


802.16 - 128kbps 


FTP File 


5 MB 


AWGN p.s.d.(AW 


-174dBm/Hz 


Pathloss exponent (y) 


3.5 


WOE ~N{[Jy dB ,Gy dB ) 


N(0dB,8dB) 


Coherent Time (Fast/Slow Fading) 


(5msec/300msec.) 


Pathloss(dB, d in meters) 


-31.5-351og I0< i + \|/ rfB 



TABLE I 

Simulation Parameters 



We will use M-LWDF-PF scheme as benchmark 
in our simulations. In this scheme at each frame the 
user maximizing the following quantity transmits. 



ai D?^(t)log(l + $P/NoW), 



(27) 



where D^ 0L (t) is the head of line packet delay and 
r,(f) is the channel capacity of user i at frame t. The 
parameter a, is a positive constant. The constant a, 
is defined as a, = — , which is referred to as 



D;" ov «,(r) 

M-LWDF-PF [7] [4]. M-LWDF-PF can be adapted 
to OFDMA systems as follows. Power is distributed 
equally to all subchannels. Starting from the first 
subchannel , the subchannel is allocated to the user 
maximizing (l27l ). Then the received rate R(t) is 
updated according to ((2]). All the subchannels are 
allocated one-by-one according to this rule. Delay 
exceeding probability is taken as 8, = 0.05 for all 
users. Delay constraint for voice, video and data 
users are 0.1, 0.4 and 1 seconds, respectively. 

Tables HT1 [TTT1 and ITVl show the effects of increasing 
the number of voice, video streaming and data 



users, respectively. In these tables we observe that 
FPQS Algorithm is better than M-LWDF algorithm 
in terms of both delay performance and data per- 
formance. With FQPS delay for voice and video 
sessions stay within acceptable bounds, while with 
M-LWDF, it exceeds the bounds for the user at 
the cell edge when V > 30 or S > 40. Besides, 
FQPS A provides at least 10 percent increase in 
total throughput. Total throughput decreases linearly 
with increasing number of voice and video users. 
Although 10 voice users adds up to 0.32 Mbps, 
adding 10 users decreases the total throughput ap- 
proximately by 1.2-1.4 Mbps. This is because voice 
has a very strict delay requirement and a voice 
session may have to be transmitted despite bad 
channel conditions. Throughput for LWDF decreases 
with a little bit slower rate but that reflects to the 
voice and video performance negatively. Log-sum 
performance of FQPSA is also better than that of 
M-LWDF, which shows that our algorithm provides 
fairness. 



V = 


10 


20 


30 


40 


50 


60 


FQPS Voice (G) 


36 


33 


30 


32 


31 


31 


FQPS Voice (B) 


88 


66 


62 


81 


82 


80 


LWDF Voice (G) 


41 


41 


41 


41 


41 


41 


LWDF Voice (B) 


89 


91 


93 


110 


133 


151 


FQPS Video (G) 


148 


132 


130 


132 


132 


131 


FQPS Video (B) 


295 


271 


250 


263 


266 


252 


LWDF Video (G) 


180 


187 


192 


200 


202 


208 


LWDF Video (B) 


384 


399 


408 


410 


423 


435 


FQPS (Mbps) 


26.9 


25.2 


23.8 


22.7 


21.4 


20.1 


FQPS log-sum 


280 


279 


278 


277 


276 


275 


LWDF (Mbps) 


23.9 


22.7 


21.6 


20.5 


19.4 


18.5 


LWDF Log- sum 


275 


278 


276 


274 


272 


269 



TABLE II 

Effects of increasing* of Voice users for 
5 = 20, D = 20. Units: Delay (msec), t.put (Mbps) 



In Table [TV] we can observe the effects of in- 
creasing the number of data users.We observe that 
delay for both voice and video streaming sessions 
stay approximately constant. Delay performance is 
much better than that of M-LWDF algorithm. Data 
performance is 10 percent better than M-LWDF. To- 
tal throughput increases with number of data users, 
but the increase diminishes as D increases. 
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s = 


10 


20 


30 


40 


50 


60 


FQPS Voice (G) 


31 


33 


39 


43 


49 


51 


FQPS Voice (B) 


57 


66 


71 


87 


93 


108 


LWDF Voice (G) 


41 


41 


41 


41 


41 


41 


LWDF Voice (B) 


86 


91 


94 


101 


142 


117 


FQPS Video (G) 


122 


132 


150 


167 


174 


190 


FQPS Video (B) 


217 


271 


275 


301 


338 


392 


LWDF Video (G) 


171 


187 


190 


194 


202 


223 


LWDF Video (B) 


395 


399 


399 


414 


463 


519 


FQPS D(Mbps) 


27.2 


25.2 


23.5 


21.5 


19.5 


17.5 


FQPS Log-sum 


281 


279 


278 


276 


274 


272 


LWDF (Mbps) 


25.1 


22.7 


20.7 


18.6 


16.7 


14.6 


LWDF log- sum 


279 


277 


277 


276 


275 


273 



TABLE III 

Effects of increasing # of Video Streaming users 
for D = 20, V = 20. Units: Delay (msec), t.put (Mbps) 



D = 


20 


30 


40 


50 


60 


70 


FQPS Voice (G) 


38 


33 


34 


34 


36 


39 


FQPS Voice (B) 


72 


66 


69 


71 


71 


91 


LWDF Voice (G) 


41 


47 


50 


57 


61 


67 


LWDF Voice (B) 


91 


106 


126 


146 


170 


182 


FQPS Video (G) 


155 


132 


143 


145 


142 


168 


FQPS Video (B) 


278 


271 


246 


270 


270 


323 


LWDF Video (G) 


187 


232 


265 


297 


335 


357 


LWDF Video (B) 


399 


460 


537 


557 


612 


656 


FQPS (Mbps) 


24.5 


25.2 


25.8 


26.1 


26.2 


26.6 


FQPS Log-sum 


146 


279 


407 


532 


654 


775 


LWDF (Mbps) 


22.7 


23.3 


23.7 


24.0 


24.1 


24.3 


LWDF Log-sum 


146 


278 


405 


530 


653 


772 



TABLE IV 

Effects of increasing # of Data users for 
5 = 20, V = 20. Units: Delay (msec), t.put (Mbps) 
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